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Abstract. We discuss exact analytical solutions of a variety of statistical models recently obtained for finite systems by a 
novel powerful mathematical method, the Laplace-Fourier transform. Among them are a constrained version of the statistical 
multifragmentation model, the Gas of Bags Model and the Hills and Dales Model of surface partition. Thus, the Laplace- 
Fourier transform allows one to study the nuclear matter equation of state, the equation of state of hadronic and quark gluon 
matter and surface partitions on the same footing. A complete analysis of the isobaric partition singularities of these models 
is done for finite systems. The developed formalism allows us, for the first time, to exactly define the finite volume analogs 
of gaseous, liquid and mixed phases of these models from the first principles of statistical mechanics and demonstrate the 
pitfalls of earlier works. The found solutions may be used for building up a new theoretical apparatus to rigorously study 
phase transitions in finite systems. The strategic directions of future research opened by these exact results are also discussed. 
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There is always a sufficient amount of facts. 
Imagination is what we lack. 

D. L Blokhintsev 

I. THEORETICAL DESCRIPTION OF PHASE TRANSITIONS IN FINITE SYSTEMS 

A rigorous theory of critical phenomena in finite systems was not built up to now. However, the experimental studies 
of phase transitions (PTs) in some systems demand the formulation of such a theory. In particular, the investigations 
of the nuclear liquid-gas PT 1 1 , 2, 3] require the development of theoretical approaches which would allow us to study 
the critical phenomena without going into the thermodynamic limit V ^ oo (V is the volume of the system) because 
such a limit does not exist due the long range Coulomb interaction. Therefore, there is a great need in the theoretical 
approaches which may shed light on the "internal mechanism" of how the PTs happen in finite systems. 

The general situation in the theory of critical phenomena for finite (small) systems is not very optimistic at the 
moment because theoretical progress in this field has been slow. It is well known that the mathematical theory of 
phase transitions was worked out by T. D. Lee and C. N. Yang fT]. Unfortunately, there is no direct generic relation 
between the physical observables and zeros of the grand canonical partition in a complex fugacity plane. Therefore, 
we know very well what are the gaseous phase and liquid at infinite volumes: mixture of fragments of all sizes and 
ocean, respectively. This is known both for pure phases and for their mixture, but, despite some limited success iQl, 
this general approach is not useful for the specific problems of critical phenomena in finite systems (see Sect. VIII 
below). 

The tremendous complexity of critical phenomena in finite systems prevented their systematic and rigorous theoret- 
ical study. For instance, even the best formulation of the statistical mechanics and thermodynamics of finite systems by 
Hill [6] is not rigorous while discussing PTs. As a result, the absence of a well established definition of the liquid and 
mixed phase for finite volumes delays the progress of several related fields, including the theoretical and experimental 
searches for the reliable signals of several PTs which are expected to exist in strongly interacting matter. Therefore, 
the task of highest priority of the theory of critical phenomena is to define the finite volume analogs of phases from 
first principles of statistical mechanics. At present it is unclear whether such definitions can be made for a general 
case, but it turns out that such finite volume definitions can be formulated for a variety of realistic nonclassical (= 
non mean-field) statistical models which are successfully used in nuclear multifragmentation and in relativistic heavy 
collsisions. 

About 25 years ago, when the theoretical foundations of nuclear multifragmentation were established, there was 
an illusion that the theoretical basis is simple and clear and, therefore, we need only the data and models which will 



describe them. The analysis of finite volume systems has proven to be very difficult. However, there was a clear way 
out of troubles by making numerical codes that are able to describe the data. This is, of course, a common way to 
handle such problems and there were many successes achieved in this way f7, . However, there is another side 

of the coin which tells us that our understanding did not change much since then. This is so because the numerical 
simulations of this level do not provide us with any proof. At best they just demonstrate something. With time the 
number of codes increased, but the common theoretical approach was not developed. This led to a bitter result - 
there are many good guesses in the nuclear multifragmentation community, but, unfortunately, little analytical work 
to back up these expectations. As a result the absence of a firm theoretical ground led to formulation of such highly 
speculative "signals" of the nuclear liquid-vapor PT as negative heat capacity L^;^, bimodality 1 10], which later on 
were disproved, in Refs fllll and fl2il . respectively. 

Thus, there is a paradoxic situation: there are many experimental data and facts, but there is no a single theoretical 
approach which is able to describe them. Similar to the searches for quark-gluon plasma (QGP) L13«l there is lack of a 
firm and rigorous theoretical approach to describe phase transitions in finite systems. 

However, our understanding of the multifragmentation phenomenon (7,^,^ was improved recently, when an exact 
analytical solution of a simplified version of the statistical multifragmentation model (SMM) 1 14, 15] was found in 
Refs. iT6lfl7ll . These analytical results not only allowed us to understand the important role of the Fisher exponent 
T on the phase structure of the nuclear liquid-gas PT and the properties of its (tri)critical point, but to calculate 
the critical indices a',l5,Y,5 of the SMM ]18] as functions of index T. The determination of the simplified SMM 
exponents allowed us to show explicitly ] 18] that, in contrast to expectations, the scaling relations for critical indices 
of the SMM differ from the corresponding relations of a well known Fisher droplet model (FDM) This exact 

analytical solution allowed us to predict a narrow range of values, 1 .799 < T < 1 .846, which, in contrast to FDM value 
Tfdm ~ 2.16, is consistent with ISiS Collaboration data ]20] and EOS Collaboration data ]21]. This finding is not only 
of a principal theoretical importance, since it allows one to find out the universality class of the nuclear liquid-gas 
phase transition, if t index can be determined from experimental mass distribution of fragments, but also it enhanced 
a great activity in extracting the value of T exponent from the data ]22]. 

It is necessary to stress that such results in principle cannot be obtained either within the widely used mean-filed 
approach or numerically. This is the reason why exactly solvable models with phase transitions play a special role 
in statistical mechanics - they are the benchmarks of our understanding of critical phenomena that occur in more 
complicated substances. They are our theoretical laboratories, where we can study the most fundamental problems of 
critical phenomena which cannot be studied elsewhere. Their great advantage compared to other methods is that they 
provide us with the information obtained directly from the first principles of statistical mechanics being unspoiled 
by mean-field or other simplifying approximations without which the analytical analysis is usually impossible. On 
the other hand an exact analytical solution gives the physical picture of PT, which cannot be obtained by numerical 
evaluation. Therefore, one can expect that an extension of the exact analytical solutions to finite systems may provide 
us with the ultimate and reliable experimental signals of the nuclear liquid-vapor PT which are established on a firm 
theoretical ground of statistical mechanics. This, however, is a very difficult general task of the critical phenomena 
theory in finite systems. 

Fortunately, we do not need to solve this very general task, but to find its solution for a specific problem of nuclear 
liquid-gas PT, which is less complicated and more realistic. In this case the straightforward way is to start from 
a few statistical models, like FDM and/or SMM, which are successful in describing the most of the experimental 
data. A systematic study of the various modifications of the FDM for finite volumes was performed by Moretto and 
collaborators fl^ and it led to a discovery of thermal reducibility of the fragment charge spectra to a determination 
of a quantitative liquid- vapor phase diagram containing the coexistence line up to critical temperature for small sy sterns 
ll24lli25l . to the generalization of the FDM for finite systems and to a formulation of the complement concept i26lE7l 
which allows one to account for finite size effects of (small) liquid drop on the properties of its vapor However, such a 
systematic analysis for the SMM was not possible until recently, when its finite volume analytical solution was found 
in(2i. 

An invention of a new powerful mathematical method I28il . the Laplace-Fourier transform, is a major theoretical 
breakthrough in the statistical mechanics of finite systems of the last decade because it allowed us to solve exactly not 
only the simplified SMM for finite volumes ]28], but also a variety of statistical surface partitions for finite clusters 
II29I1 and to find out their surface entropy and to shed light on a source of the Fisher exponent T. It was shown 128,1 
that for finite volumes the analysis of the grand canonical partition (GCP) of the simplified SMM is reduced to the 
analysis of the simple poles of the corresponding isobaric partition, obtained as a Laplace-Fourier transform of the 
GCP. Such a representation of the GCP allows one not only to show from first principles that for finite systems 
there exist the complex values of the effective chemical potential, but to define the finite volume analogs of phases 



straightforwardly. Moreover, this method allows one to include into consideration all complicated features of the 
interaction (including the Coulomb one) which have been neglected in the simplified SMM because it was originally 
formulated for infinite nuclear matter. Consequently, the Laplace-Fourier transform method opens a principally new 
possibility to study the nuclear liquid-gas phase transition directly from the partition of finite system without taking 
its thermodynamic limit. Now this method is also applied |30] to the finite volume formulation of the Gas of Bags 
Model (GBM) 1 3 1 ] which is used to describe the PT between the hadronic matter and QGP. Thus, the Laplace-Fourier 
transform method not only gives an analytical solution for a variety of statistical models with PTs in finite volumes, 
but provides us with a common framework for several critical phenomena in strongly interacting matter. Therefore, 
it turns out that further applications and developments of this method are very promising and important not only for 
the nuclear multifragmentation community, but for several communities studying PTs in finite systems because this 
method may provide them with the firm theoretical foundations and a common theoretical language. 

It is necessary to remember that further progress of this approach and its extension to other communities cannot be 
successfully achieved without new theoretical ideas about formalism it-self and its applications to the data measured in 
low and high energy nuclear collisions. Both of these require essential and coherent efforts of two or three theoretical 
groups working on the theory of PTs in finite systems, which, according to our best knowledge, do not exist at the 
moment either in multifragmentation community or elsewhere. Therefore, the second task of highest priority is to 
attract young and promising theoretical students to these theoretical problems and create the necessary manpower to 
solve the up coming problems. Otherwise the negative consequences of a complete dominance of experimental groups 
and numerical codes will never be overcome and a good chance to build up a common theoretical apparatus for a few 
PTs will be lost forever. If this will be the case, then an essential part of the nuclear physics associated with nuclear 
multifragmentation will have no chance to survive in the next years. 

Therefore, the first necessary step to resolve these two tasks of highest priority is to formulate the up to day achieve- 
ments of the exactly solvable models and to discuss the strategy for their further developments and improvements 
along with their possible impact on transport and hydrodynamic approaches. For these reasons the paper is organized 
as follows: in Sect. II we formulate the simplified SMM and present its analytical solution in thermodynamic limit; in 
Sect. Ill we discuss the necessary conditions for PT of given order and their relation to the singularities of the isobaric 
partition and apply these findings to the simplified SMM; Sect. IV is devoted to the SMM critical indices as the func- 
tions of Fisher exponent t and their scaling relations; the Laplace-Fourier transform method is presented in Sect. V 
along with an exact analytical solution of the simplified SMM which has a constraint on the size of largest fragment, 
whereas the analysis of its isobaric partition singularities and the meaning of the complex values of free energy are 
given in Sect. VI; Sect. VII and VIII are devoted to the discussion of the case without PT and with it, respectively; at 
the end of Sect. VIII there is a discussion of the Chomaz and Gulminelli's approach to bimodality |5|; in Sect. IX we 
discuss the finite volume modifications of the Gas of Bags, i.e. the statistical model describing the PT between hadrons 
and QGP, whereas in Sect. X we formulate the Hills and Dales Model for the surface partition and present the limit of 
the vanishing amplitudes of deformations; and, finally, in Sect. XI we discuss the strategy of future research which is 
necessary to build up a truly microscopic kinetics of phase transitions in finite systems. 

II. STATISTICAL MULTIFRAGMENTATION IN THERMODYNAMIC LIMIT 

The system states in the SMM are specified by the multiplicity sets {n^} (n^. = 0, 1,2, ...) of A:-nucleon fragments. The 
partition function of a single fragment with k nucleons is (ij]: y0^.(r) = V {mTk/ln)^^^ Zk , where A: = 1,2, ...,A (A is 
the total number of nucleons in the system), V and T are, respectively, the volume and the temperature of the system, m 
is the nucleon mass. The first two factors on the right hand side (r.h.s.) of the single fragment partition originate from 
the non-relativistic thermal motion and the last factor, Zk, represents the intrinsic partition function of the A:-nucleon 
fragment. Therefore, the function ^k{T) is a phase space density of the k-nucleon fragment. For k—\ (nucleon) we 
take z\ = 4 (4 internal spin-isospin states) and for fragments with > 1 we use the expression motivated by the liquid 
drop model (see details in Ref. Hi]): Zk — ^wi^fk/T), with fragment free energy 

fk = -W{T) k+a{T) ^2/3 + (T + 3/2)rin^ , (1) 

with W{T) = Wo + T^/Eo. Here Wo = 16 MeV is the bulk binding energy per nucleon. T^/Eq is the contribution of the 
excited states taken in the Fermi-gas approximation (Eq = 16 MeV). Ct{T) is the temperature dependent surface tension 
parameterized in the following relation: a{T) = a{T)\sMM = OoliT^ - T^)/{T^ + r2)]V4^ with Go = 18 MeV and 
= 18 MeV (a — at T > Tc). The last contribution in Eq. Q involves the famous Fisher's term with dimensionless 
parameter T. 



It is necessary to stress that the SMM parametrization of the surface tension coefficient is not a unique one. For 
instance, the FDM successfully employs another one g{T)\pum = <7o[l — T/Tc]. As we shall see in Sect. IV the 
temperature dependence of the surface tension coefficient in the vicinity of the critical point will define the critical 
indices of the model, but the following mathematical analysis of the SMM is general and is valid for an arbitrary C7(r) 
function. 

The canonical partition function (CPF) of nuclear fragments in the SMM has the following form: 
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In Eq. ^ the nuclear fragments are treated as point-like objects. However, these fragments have non-zero proper 
volumes and they should not overlap in the coordinate space. In the excluded volume (Van der Waals) approximation 
this is achieved by substituting the total volume V in Eq. (|3 by the free (available) volume Vf = V --bY^kkrik, 
where b = l/po (Po = 0. 16 fm^^ is the normal nuclear density). Therefore, the corrected CPF becomes: Za {V, T) = 
Z'^iy — bA, T). The SMM defined by Eq. Q was studied numerically in Refs. lITillTsll . This is a simplified version of 
the SMM, since the symmetry and Coulomb contributions are neglected. However, its investigation appears to be of 
principal importance for studies of the nuclear liquid-gas phase transition. 

The calculation of ZA(y, T) is difficult due to the constraint Y.k^'^k = A. This difficulty can be partly avoided by 
evaluating the grand canonical partition (GCP) 



^{VJ,^) = £exp(^^)zA(y,r)0(V-M), 



(3) 



A=0 



where jj. denotes a chemical potential. The calculation of 3f is still rather difficult. The summation over {n^ } sets in Za 
cannot be performed analytically because of additional A-dependence in the free volume Vf and the restriction Vf > 0. 
The presence of the theta-function in the GCP guarantees that only configurations with positive value of the free 
volume are counted. However, similarly to the delta function restriction in Eq. 0, it makes again the calculation of 
iF(y, T,IJ.) (13 to be rather difficult. This problem was resolved 1 16. il7ll by p erforming the Laplace transformation of 
iF(y, T,n). This introduces the so-called isobaric partition function (IP) f3V\: 
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After changing the integration variable V V\ the constraint of ©-function has disappeared. Then all rik were summed 
independently leading to the exponential function. Now the integration over V' in Eq. can be done resulting in 
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Here we have introduced the shifted chemical potential ji = ji + W{T). From the definition of pressure in the grand 
canonical ensemble it follows that, in the thermodynamic limit, the GCP of the system behaves as 
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An exponentially over V increasing part of if°(y, in the right-hand side of Eq. Q generates the rightmost 
singularity s* of the function ^{s,T,jj.), because for s < p{T,jj.)/T the y-integral for 3f{s,T,iJ.) diverges at its 
upper limit. Therefore, in the thermodynamic limit, V ^ °° the system pressure is defined by this rightmost singularity, 
s*{T,ii),of IP S'is,!, 11) (Hi: 

p{T,ii) ^Ts*{T,ii). (8) 

Note that this simple connection of the rightmost i-singularity of iF, Eq. 0, to the asymptotic, V ^ °°, behavior of 
jF, Eq. Q, is a general mathematical property of the Laplace transform. Due to this property the study of the system 
behavior in the thermodynamic limit V ^ °° can be reduced to the investigation of the singularities of iF. 



III. SINGULARITIES OF ISOBARIC PARTITION AND PHASE TRANSITIONS 



The IP, Eq. has two types of singularities: 1) the simple pole singularity defined by the equation 



(9) 



2) the singularity of the function ^{s,T,jj.) it-self at the point si where the coefficient in linear over k terms in the 
exponent is equal to zero, 



The simple pole singularity corresponds to the gaseous phase where pressure is determined by the equation 
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The singularity i/(r,/i) of the function ,^{s, T,ii) (|6j defines the liquid pressure 



(10) 
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In the considered model the liquid phase is represented by an infinite fragment, i.e. it corresponds to the macroscopic 
population of the single mode k = °o. Here one can see the analogy with the Bose condensation where the macroscopic 
population of a single mode occurs in the momentum space. 

In the -regions where p. < bpg{T,iJ.) the gas phase dominates (pg > pi), while the liquid phase corresponds 
to /I > bpg{T^fx). The liquid-gas phase transition occurs when two singularities coincide, i.e. Sg{T,ii) — A 
schematic view of singular points is shown in Fig. la for T < T^, i.e. when a > 0. The two-phase coexistence region 
is therefore defined by the equation 



Pi{T,ii) = Pg{T,fi) , i.e., A = bpg{T,ii) . 



(13) 



One can easily see that ,^{s,T,pL) is monotonously decreasing function of s. The necessary condition for the phase 
transition is that this function remains finite in its singular point si = fi/Tb: 



(14) 



The convergence of is determined by T and d. At T = the condition (I14> requires g{T) > 0. Otherwise, 
■^{si,T,ji) = oo and the simple pole singularity Sg{T,jj.) (|9} is always the rightmost ^-singularity of (0} (see 
Fig. lb). At r > Tc, where g{T)\smm = 0, the considered system can exist only in the one-phase state. It will be 
shown below that for t > 1 the condition (tT4ll can be satisfied even at (j{T) = 0. 

At r < Tc the system undergoes the 1-st order phase transition across the line jj.* — IJ.*{T) defined by Ea.( ll3> . Its 
explicit form is given by the expression: 
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The points on the line li*{T) correspond to the mixed phase states. First we consider the case T = —1.5 because it is 
the standard SMM choice. 

The baryonic density is found as {dp/dii)T and is given by the following formulae in the liquid and gas phases 



Pi 



djj. 



respectively. Here the function p,v/ is defined as 



■r\ 3/2 
ml ^ ' 

'in 



Z\ exp 



1 
b 



l^-bp^ 



dpi 



Pid 



^+bpid 



k=2 



1-T 



exp I 



' {jl-bpg)k-ak^/^^ 



(16) 



(17) 



s 

a) b) 

FIGURE 1. Schematic view of singular points of the Isobaric Partition, Eq. Js}, at T < (a) and T > Tc (b). Full lines show 
,!^{s,T,ll) as a function of at fixed T and /i, /Xi < ^2 < M3 < M4- Dots and asterisks indicate the simple poles (Sg) and the 
singularity of function JP' it-self (s/). At fiT, = p.*{T) the two singular points coincide signaling a phase transition. 
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FIGURE 2. Phase diagram in T — p plane for T = —1.5 (left), T = 1.1 (middle) and T = 2.1 (right). The mixed phase is 
represented by the extended region. Liquid phase (shown by crosses) exists at density p = Po. Point C in the left panel is the critical 
point, whereas in the middle panel it is the tricritical point. For T > 2 (right panel) the PT exists for all temperatures T >0. 



Due to the condition ( I13> this expression is simplified in the mixed phase: 
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This formula clearly shows that the bulk (free) energy acts in favor of the composite fragments, but the surface term 
favors single nucleons. 

Since at a > the sum in Eq. ( I18> converges at any T, p,y is finite and according to Eq. (I16t Pg < 1/b. Therefore, 
the baryonic density has a discontinuity Ap = p/ — Pg > across the line u*(T) O for any T. The discontinuities 
take place also for the energy and entropy densities. The phase diagram of the system in the (r,p)-plane is shown 
in the left panel of Fig. 2. The line l-L*{T) il5i corresponding to the mixed phase states is transformed into the finite 
region in the (r,p)-plane. As usual, in this mixed phase region of the phase diagram the baryonic density p and the 
energy density are superpositions of the corresponding densities of liquid and gas: 



p = X Pi + (1 - X) Pg , e = A £/ + (1 - A) 



(19) 



Here A(0<A<l)isa fraction of the system volume occupied by the liquid inside the mixed phase, and the partial 
energy densities for (/ = l,g) can be found from the thermodynamic identity flfill : 

Inside the mixed phase at constant density p the parameter A has a specific temperature dependence shown in the 
left panel of Fig. 3: from an approximately constant value p/po at small T the function A {T) drops to zero in a narrow 
vicinity of the boundary separating the mixed phase and the pure gaseous phase. This corresponds to a fast change 
of the configurations from the state which is dominated by one infinite liquid fragment to the gaseous multifragment 
configurations. It happens inside the mixed phase without discontinuities in the thermodynamical functions. 

An abrupt decrease of X{T) near this boundary causes a strong increase of the energy density as a function of 
temperature. This is evident from the middle panel of Fig. 3 which shows the caloric curves at different baryonic 
densities. One can clearly see a leveling of temperature at energies per nucleon between 10 and 20 MeV. As a 
consequence this leads to a sharp peak in the specific heat per nucleon at constant density, Cp{T) = {de/dT)p /p , 
presented in Fig. 3. A finite discontinuity of Cp{T) arises at the boundary between the mixed phase and the gaseous 
phase. This finite discontinuity is caused by the fact that A [T) = 0, but {dX /dT)p ^ at this boundary (see Fig. 3). 

It should be emphasized that the energy density is continuous at the boundary of the mixed phase and the gaseous 
phase, hence the sharpness of the peak in Cp is entirely due to the strong temperature dependence of A(r) near this 
boundary. Moreover, at any p < the maximum value of Cp remains finite and the peak width in Cp{T)js nonzero 
in the thermodynamic limit considered in our study. This is in contradiction with the expectation of Refs. |l4lfl5| that 
an infinite peak of zero width will appear in Cp{T) in this limit. Also a comment about the so-called "boiling point" 
is appropriate here. This is a discontinuity in the energy density e(r) or, equivalently, a plateau in the temperature as 
a function of the excitation energy. Our analysis shows that this type of behavior indeed happens at constant pressure, 
but not at constant density! This is similar to the usual picture of a liquid-gas phase transition. In Refs. lll4iri5ll a rapid 
increase of the energy density as a function of temperature at fixed p near the boundary of the mixed and gaseous 
phases (see the right panel of Fig. 3) was misinterpreted as a manifestation of the "boiling point". 

New possibilities appear at non-zero values of the parameter T. At < T < 1 the qualitative picture remains the 
same as discussed above, although there are some quantitative changes. For T > 1 the condition ( I14> is also satisfied at 
T > Tc where o{T)\smm — 0. Therefore, the liquid-gas phase transition extends now to all temperatures. Its properties 
are, however, different for T > 2 and for T < 2 (see Fig. 2). If T > 2 the gas density is always lower than 1 /b as p,y 
is finite. Therefore, the liquid-gas transition at T > 71. remains the 1-st order phase transition with discontinuities of 
baryonic density, entropy and energy densities (right panel in Fig. 2) . 



IV. THE CRITICAL INDICES AND SCALING RELATIONS OF THE SMM 

The above results allow one to find the critical exponents a', )3 and / of the simplified SMM. These exponents describe 
the temperature dependence of the system near the critical point on the coexistence curve ji* — pi*{T) (I13> . where the 
effective chemical potential vanishes v = pi*{T)+W{T) -bp{T,ii*{T)) ^0 
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where Ap = p/ — Pg defines the order parameter, Cp = ^ ^ denotes the specific heat at fixed particle density and 

Kt = ^ (^) isothermal compressibility. The shape of the critical isotherm for p < p^. is given by the critical 

index 5 (the tilde indicates e = hereafter) 
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The calculation of a and a' requires the specific heat Cp. With the formula ll34ll 
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FIGURE 3. Left panel: Volume fraction A(r) of the liquid inside the mixed phase is shown as a function of temperature for 
fixed nucleon densities p /po = 1/6, 1/3, 1/2, 2/3, 5/6 (from bottom to top) and t = —1.5. 

Middle panel: Temperature as a function of energy density per nucleon (caloric curve) is shown for fixed nucleon densities 
p/Po — 1/6, 1/3, 1/2,2/3 and T = —1.5. Note that the shape of the model caloric curves is very similar to the experimental 
finding 1 32], although our estimates for the excitation energy is somewhat larger due to oversimplified interaction. For a quantitative 
comparison between the simplified SMM the full SMM interaction should be accounted for. 

Right panel: Specific heat per nucleon as a function of temperature at fixed nucleon density p/po = 1/3. The dashed line shows 
the finite discontinuity of Cp (r) at the boundary of the mixed and gaseous phases for T = — 1 .5. 



one obtains the specific heat on the PT curve by replacing the partial derivatives by the total ones f33\. The latter 
can be done for every state inside or on the boundary of the mixed phase region. For the chemical potential /X* (T) = 



bp* {T)-W{T) one gets 



dV(i) , d%(n 



AT- 



Here the asterisk indicates the condensation line (v — 0) 



hereafter. Fixing p — pi = l/b one finds c*p^{T) = T ' and, hence, obtains a = a' = 0. To calculate j3, / 
and 5 the behavior of the series 

Eq(e,v) = £ k'^-'eT/ (26) 

k=2 

should be analyzed for small positive values of £ and — v (A = Oo/Tc)- In the limit e — ^ the function Eq(e,0) remains 
finite, if T > ^ + 1, and diverges otherwise. For T — q+l this divergence is logarithmic. The case T < ^ + 1 is analyzed 
in some details, since even in Fisher's papers it was performed incorrectly. 
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one can prove f 181 that in the limit e ^ the series on the r h. s. of i26i 
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The assumption q — t > — 1 is sufficient to guarantee the convergence of the integral at its lower limit. Using this 
representation, one finds the following general results fisll 



Eq(e,0) 



eii'^ 1 1). ifT<^+l. 
ln|£|, ifT = <7+l, and Sq(0,v) 



,0 



if T > ^ + 1 , 



yT-q-l^ ifT<^+l, 
In I V I , if T = ^ + 1 , 

if T > ^ + 1 , 



(28) 



which allowed us to find out the critical indices of the SMM (see Table 1). 
In the special case ^ = 2a the well-known exponent inequalities proven for real gases by 



Fisher[36] : a' + ljS + Y > 2, 



(29) 



TABLE 1. Critical exponents of the SMM and FDM as functions of Fisher index T for the 
general parametrization of the surface energy (y(T)0 e^k" withe = (Tc-T)/Tc. 
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Griffiths [37] : a' + j3(l+5) >2, (30) 
LibermanpS]: / + j3(l-5) >0, (31) 

are fulfilled exactly for any T. (The corresponding exponent inequalities for magnetic systems are often called Rush- 
brooke's, Griffiths' and Widom's inequalities, respectively.) For ^ > 2a, Fisher's and Griffiths' exponent inequalities 
are fulfilled as inequalities and for ^ < 2a they are not fulfilled. The contradiction to Fisher's and Griffiths' exponent 
inequalities in this last case is not surprising. This is due to the fact that in the present version of the SMM the critical 
isochore p = Pc = pi lies on the boundary of the mixed phase to the liquid. Therefore, in expression (2.13) in Ref. 
lis^ for the specific heat only the liquid phase contributes and, therefore. Fisher's proof of Ref. ll3^ following (2.13) 
cannot be applied for the SMM. Thus, the exponent inequalities (I29> and i30\ have to be modified for the SMM. Using 
results of Tabled one finds the following scaling relations 

a' + 2l5 + Y = - and a' + j3(l + 5) = - . (32) 

(7 (7 

Liberman's exponent inequality ( 13 11 1 is fulfilled exactly for any choice of ^ and a. 

Since the coexistence curve of the SMM is not symmetric with respect to p = Pc, it is interesting with regard to 
the specific heat to consider the difference Acp{T) = c*p^{T) ~c*p^{T), following the suggestion of Ref. [35J . Using 
Eq. for gas and liquid and noting that 1 / p* ~b — 1/p,*;, one obtains a specially defined index a'^ from the most 
divergent term for ^ > 1 



i 

1, I J ) — -, — r ^ — u,.. — s ^ 

" P,-,(r) dr- ' \ 2-|(o + 2-t), ifi>T + i, 



^An-J^'-^ ^ , "^<"'+'- ,33) 



Then it is a'^ > for ^/a <2. Thus, approaching the critical point along any isochore within the mixed phase region 
except for p = p^. = 1//? the specific heat diverges for i^/ a < 2 as defined by a'^ and remains finite for the isochore 
p = Pc — l/b. This demonstrates the exceptional character of the critical isochore in this model. 

In the special case that ^ = I one finds a( =2 — 1/(7 for T < 1 + 2a and a'^ ~ — )3 for T > 1+2(7. Therefore, 
using ttj instead of a', the exponent inequalities ( I29> and ( I30t are fulfilled exactly if ^ > 1 and T<(7+lorifi^ = l 
and T < 2(7 + 1. In all other cases ( I29l l and ( I30t are fulfilled as inequalities. Moreover, it can be shown that the SMM 
belongs to the universality class of real gases for C > 1 and f > a+l. 

The comparison of the above derived formulae for the critical exponents of the SMM for C = 1 with those obtained 
within the FDM (Eqs. 51-56 in L191 ) shows that these models belong to different universality classes (except for the 
singular case t — 2). 

Furthermore, one has to note that for^ = l, a <l/2 and 1 + a < T < 1 + 2(7 the critical exponents of the SMM 
coincide with those of the exactly solved one-dimensional FDM with non-zero droplet-volumes 1 35]. 
For the usual parameterization of the SMM (Jl one obtains with 1^ — 5/4- and a = 2/3 the exponents 

/ ^' « 15 ,^ , , 15 / 3\ . T-1 

«v = i „ , . ^ = ^ 2-T , y = It- \ 5 = ——. (34) 

[ ^T-3, if T> I 8 4 V 27 2-T 



Thus, Fisher suggestion to use a'^ instead of a' allows one to "save " the exponential inequalities, however, it is not a 
final solution of the problem. 



The critical indices of the nuclear Uquid-gas PT were determined from the multifragmentation of gold nuclei |3S 
and found to be close to those ones of real gases. The method used to extract the critical exponents j3 and / in Ref. 
II39'] was, however, found to have large uncertainties of about 25 per cents f4^. Nevertheless, those results allow us to 
estimate the value of T from the experimental values of the critical exponents of real gases taken with large error bars. 
Using the above results we generalized llisll the exponent relations of Ref. l35ll 

for arbitrary a and ^. Then, one obtains with l41ll j3 = 0.32 — 0.39, / = 1.3 — 1.4 and 5=4 — 5 the estimate 
T — 1.799 — 1.846. This demonstrates also that the value of T is rather insensitive to the special choice of j3 , / 
and 5, which leads to a!. = 0.373 — 0.461 for the SMM. Theoretical values for /3 , / and 8 for Ising-like systems 
within the renormalized 0'* theory |42] lead to the narrow range T = 1.828 ±0.001 . The values of j3 , / and 5 indices 
for nuclear matter and percolation of two- and three-dimensional clusters are reviewed in 123.1 . 

There was a decent try to study the critical indices of the SMM numerically ll44ll . The version V2 of Ref. ll44ll 
corresponds precisely to our model with T = 0, ^=5/4 and (7 = 2/3, but their results contradict to our analysis. Their 
results for version V3 of Ref. L44] are in contradiction with our proof presented in Ref. L161 . There it was shown that 
for non-vanishing surface energy (as in version V3) the critical point does not exist at all. The latter was found in 
f47] for the finite system and the critical indices were analyzed. Such a strange result, on one hand, indicates that the 
numerical methods used in Ref. |44] are not self-consistent, and, on the other hand, it shows an indispensable value of 
the analytical calculations, which can be used as a test problem for numerical algorithms. 

It is widely believed that the effective value of T defined as Teff = — (9 In ( e ) /(9 In ^ attains its minimum at the critical 
point (see references in L21.1 ). This has been shown for the version of the FDM with the constraint of sufficiently small 
surface tension a = for T >Tc ll43ll and also can be seen easily for the SMM. Taking the SMM fragment distribution 
ni,{e) = g{T)k-^ exp[^k - ^k"] ~ k-'^>=« one finds 

V (Tfl(£) 

Teff = T:--k+ -^j^k" ^ T = min(Teff ) , (36) 

where the last step follows from the fact that the inequalities a(e) > , V < become equalities at the critical point 
V = fl(0) = 0. Therefore, the SMM predicts that the minimal value of Teff coiTesponds to the critical point. 

In the E900 k^+Pm multifragmentation experiment i2(il the ISiS collaboration measured the dependence of 
Teff upon the excitation energy and found the minimum value min(Teff) = 1.9 (Fig. 5 of Ref. |20] ). Also the EOS 
collaboration |21] performed an analysis of the minimum of Teff on AuH-C multifragmentation data. The fitted Teff, 
plotted in Fig. 1 l.b of Ref. |21] versus the fragment multiplicity, exhibits a minimum in the range min(Tejf) = 1.8 — 1.9 . 
Both results contradict the original FDM 1 19], but agree well with the above estimate of T for real gases and for Ising- 
like systems in general. 



V. CONSTRAINED SMM IN FINITE VOLUMES 

Despite the great success, the application of the exact solution 

EES El 

to the description of experimental data is 
limited because this solution corresponds to an infinite system and due to that it cannot account for a more complicated 
interaction between nuclear fragments. Therefore, it was necessary to extend the exact solution 1 16, 17, 18] to finite 
volumes. It is clear that for the finite volume extension it is necessary to account for the finite size and geometrical 
shape of the largest fragments, when they are comparable with the system volume. For this one has to abandon the 
arbitrary size of largest fragment and consider the constrained SMM (CSMM) in which the largest fragment size is 

K(V) 

explicitly related to the volume V of the system. Thus, the CSMM assumes a more strict constraint ^ k nj^ — A , 

k 

where the size of the largest fragment K{V) — aV /b cannot exceed the total volume of the system (the parameter 
a < 1 is introduced for convenience). The case of the fixed size of the largest fragment, i.e. KiV) — Const, analyzed 
numerically in Ref. i45ll is also included in our treatment. A similar restriction should be also applied to the upper 
limit of the product in all partitions Z^^iV.T), Za{V,T) and S'{V,T,ji) introduced above (how to deal with the real 
values of K{V), see later). Then the model with this constraint, the CSMM, cannot be solved by the Laplace transform 
method, because the volume integrals cannot be evaluated due to a comphcated functional V-dependence. However, 



the CSMM can be solved analytically with the help of the following identity |2l 



(37) 



which is based on the Fourier representation of the Dirac 5-function. The representation ( I37t allows us to decouple 
the additional volume dependence and reduce it to the exponential one, which can be dealt by the usual Laplace 
transformation in the following sequence of steps 



ir(A,r,jU) EE £ dV e-'^^ ^{V,T,fl)^ j^dV' j d£, J 
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After changing the integration variable V ^V' — b Y, kuj^, the constraint of ©-function has disappeared. Then 

k 

all iiif were summed independently leading to the exponential function. Now the integration over V' in Eq. (I38> can be 
straightforwardly done resulting in 



d^ 



ir(A,r,M) 



where the function A) is defined as follows 

[tl-lbT)k 

k=\ 



dr\ 

2% X-iri - .^{^,X-iri) 



mT 
2n 



k=2 



(ll+W-lhT)k-ak^l^ 

T 



(39) 



(40) 



As usual, in order to find the GCP by the inverse Laplace transformation, it is necessary to study the structure of 
singularities of the isobaric partition (I40> . 

VI. ISOBARIC PARTITION SINGULARITIES AT FINITE VOLUMES 

The isobaric partition (I40t of the CSMM is, of course, more complicated than its SMM analog f\3, because for 
finite volumes the structure of singularities in the CSMM is much richer than in the SMM, and they match in the limit 
V ^ oa only. To see this let us first make the inverse Laplace transform: 



'di M_lf^_ 

J iK J Ini X - /Tj 



X-'" 



y 



(41) 



noindent where the contour A -integral is reduced to the sum over the residues of all singular points A = A„ + ; 77 with 
n = 1,2,.., since this contour in the complex A-plane obeys the inequality x > max(/?e{A„}). Now both remaining 
integrations in ( 14 U can be done, and the GCP becomes 



ir(y,r,M) 



.A„y 



-1 



(42) 



i.e. the double integral in (I41> simply reduces to the substitution E, ^ V in the sum over singularities. This is a 
remarkable result which was formulated in Ref. |28] as the following theorem: if the Laplace-Fourier image of the 
excluded volume GCP exists, then for any additional V -dependence of ^{V, A„) or <j)k{T) the GCP can be identically 
represented by Eq. \42\ . 

The simple poles in (I41> are defined by the equation 



A„ - jr(y,A„) 



(43) 



In contrast to the usual SMM JH [H the singularities A.„ are (i) are volume dependent functions, if K{V) is not 
constant, and (ii) they can have a non-zero imaginary part, but in this case there exist pairs of complex conjugate roots 
of j43> because the GCP is real. 

Introducing the real R„ and imaginary /„ parts of A„ — R„ + il„, we can rewrite Eq. ( l43t as a system of coupled 
transcendental equations 

K{V) ^ 

Rn = L ^kiT) cos(/„M) , (44) 

k=\ 

- - E ~^k(J) ^m[lnhk) , (45) 

k=\ 

where we have introduced the set of the effective chemical potentials V„ = v(A„) with v(A) = [I +W(r) —XbT, 

and the reduced distributions = (f|^)'ziexp(-W^(r)/r) and ^k>\{T) = {^Y k-'' exp(-c7(r) k^l^ /T) for 

convenience. 

Consider the real root {Rq > 0,/o = 0), first. For /„ = /q = the real root Rq exists for any T and jU. Comparing 
Ro with the expression for vapor pressure of the analytical SMM solution [16. 17.1 shows that TRq is a constrained 
grand canonical pressure of the gas. As usual, for finite volumes the total mechanical pressure 0, E^l differs from 
TRq. Equation ( I45t shows that for /„>o 7^ the inequality cos{I„bk) < 1 never become the equality for all A;- values 
simultaneously. Then from Eq. (I44t one obtains (« > 0) 

^« < L '^^(^) ^ Rn<RQ, (46) 

<r=l 

where the second inequality (I46> immediately follows from the first one. In other words, the gas singularity is always 
the rightmost one. This fact plays a decisive role in the thermodynamic limit V ^ °°. 

The interpretation of the complex roots A„>o is less straightforward. According to Eq. (l42ll . the GCP is a superpo- 
sition of the states of different free energies —X„VT. (Strictly speaking, —X„VT has a meaning of the change of free 
energy, but we will use the traditional term for it.) For n > the free energies are complex. Therefore, —X„yQT is the 
density of free energy. The real part of the free energy density, —R„T, defines the significance of the state's contri- 
bution to the partition: due to ( I46> the largest contribution always comes from the gaseous state and has the smallest 
real part of free energy density. As usual, the states which do not have the smallest value of the (real part of) free 
energy, i. e. —R„yQT, are thermodynamically metastable. For infinite volume they should not contribute unless they 
are infinitesimally close to —RqT, but for finite volumes their contribution to the GCP may be important. 

As one sees from (I44> and j45> . the states of different free energies have different values of the effective chemical 
potential v„, which is not the case for infinite volume where there exists a single value for the effective 

chemical potential. Thus, for finite V the states which contribute to the GCP ( I42> are not in a true chemical equilibrium. 

The meaning of the imaginary part of the free energy density becomes clear from (I44> and (145 > |46]: as one can see 
from (|43 the imaginary part /„>o effectively changes the number of degrees of freedom of each fc-nucleon fragment 
(k < K{V)) contribution to the free energy density —R„^oT- It is clear, that the change of the effective number of 
degrees of freedom can occur virtually only and, if A„>o state is accompanied by some kind of equilibration process. 
Both of these statements become clear, if we recall that the statistical operator in statistical mechanics and the quantum 
mechanical convolution operator are related by the Wick rotation |47|. In other words, the inverse temperature can 
be considered as an imaginary time. Therefore, depending on the sign, the quantity I^bT = \ ^ that appears in the 
trigonometric functions of the equations ( I44> and ( I45t in front of the imaginary time \/T can be regarded as the 
inverse decay/formation time T„ of the metastable state which corresponds to the pole A,„>o (for more details see next 
sections). 

This interpretation of T„ naturally explains the thermodynamic metastability of all states except the gaseous one: 
the metastable states can exist in the system only virtually because of their finite decay/formation time, whereas the 
gaseous state is stable because it has an infinite decay/formation time. 



VII. NO PHASE TRANSITION CASE 



It is instructive to treat the effective chemical potential v{X) as an independent variable instead of jU. In contrast to 
the infinite V , where the upper limit v < defines the liquid phase singularity of the isobaric partition and gives 




0.0 0.2 0.4 0.6 0.8 1.0 5 10 IS 20 25 

Ijb T (MeV) 

FIGURE 4. Left panel: A graphical solution of Eq. {45} for T = 10 MeV and T = 1.825. The l.h.s. (straight line) and r.h.s. of 
Eq. 145 > (all dashed curves) are shown as the function of dimensionless pa ram eter /i b for the three values of the largest fragment 
size K(V). The intersection point at (0; 0) corresponds to a real root of Eq. <43> . Each tangent point with the straight line generates 
two complex roots of <43> . 

Right panel: Each curve separates the T — Re{v„) region of one real root of Eq. <43t (below the curve), three complex roots (at the 
curve) and four and more roots (above the curve) for three values of K{V) and the same parameters as in the left panel. 



the pressure of a liquid phase pi{T,ii) = TRo\v^oo = (/i + W{T))/b iH^IitIi . for finite volumes and finite K{V) the 
effective chemical potential can be complex (with either sign for its real part) and its value defines the number and 
position of the imaginary roots {A„>o} in the complex plane. Positive and negative values of the effective chemical 
potential for finite systems were considered |33| within the Fisher droplet model, but, to our knowledge, its complex 
values have never been discussed. From the definition of the effective chemical potential v(A) it is evident that its 
complex values for finite systems exist only because of the excluded volume interaction, which is not taken into 
account in the Fisher droplet model 1 19]. However, a recent study of clusters of the d = 2 Ising model within the 
framework of FDM (see the corresponding section in Ref. [23J) shows that the excluded volume correction improves 
essentially the description of the thermodynamic functions. Therefore, the next step is to consider the complex values 
of the effective chemical potential and free energy for the excluded volume correction of the Ising clusters on finite 
lattices. 

As it is seen from Fig. 1, the rh.s. of Eq. j45> is the amplitude and frequency modulated sine-like function of 
dimensionless parameter /„ b. Therefore, depending on T and Re{v) values, there may exist no complex roots {A„>o}, 
a finite number of them, or an infinite number of them. In Fig. 1 we showed a special case which corresponds to exactly 
three roots of Eq. (I43> for each value of K(y): the real root (/q = 0) and two complex conjugate roots (±/i). Since 
the r.h.s. of ( l45l is monotonously increasing function of Re{v), when the former is positive, it is possible to map the 
T —Re{v) plane into regions of a fixed number of roots of Eq. ( I43t . Each curve in Fig. 2 divides the T — Re{v) plane 
into three parts: for 7?e(v)-values below the curve there is only one real root (gaseous phase), for points on the curve 
there exist three roots, and above the curve there are four or more roots of Eq. ( I43> . 

For constant values of K{V) = K the number of terms in the r.h.s. of J45> does not depend on the volume and, 
consequently, in thermodynamic limit V ^ °° only the rightmost simple pole in the complex -plane survives out of a 
finite number of simple poles. According to the inequality ( I46t . the real root Aq is the rightmost singularity of isobaric 
partition (13 9> . However, there is a possibility that the real parts of other roots A„>o become infinitesimally close to Rq, 
when there is an infinite number of terms which contribute to the GCP \A2l . 

Let us show now that even for an infinite number of simple poles in (I42> only the real root Xq survives in the limit 
y ^ oo. For this purpose consider the limit Re{Vn) ^ T . In this limit the distance between the imaginary parts of the 
nearest roots remains finite even for infinite volume. Indeed, for Re{v„) ^ T the leading contribution to the r.h.s. of 
(145 > corresponds to the harmonic with k = K, and, consequently, an exponentially large amplitude of this term can be 
only compensated by a vanishing value of sin {I„bK), i.e. I„bK = 7rn + d„ with |5„| ^ ;r (hereafter we will analyze 
only the branch /„ > 0), and, therefore, the corresponding decay /formation time t„ w K[7:nT] ^ ^ is volume independent. 



Keeping the leading term on the r.h.s. of J45t and solving for 5„, one finds 



5„, with 5„ 



Re(v„)K 
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(47) 
(48) 



where in the last step we used Eq. (I44> and condition |5„| ^ TT. Since for V ^ oo all negative values of R„ cannot 
contribute to the GCP (I42> . it is sufficient to analyze even values of n which, according to ( I48> . generate R„ > 0. 

Since the inequality (I46> can not be broken, a single possibility, when A„>o pole can contribute to the partition ( I42> . 
corresponds to the case R„ Rq—Q^ for some finite n. Assuming this, we find Re{v{Xn)) 7?e(v(Ao)) for the same 
value of ji. Substituting these results into equation J44> . one gets 



k=\ 



nnk 
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(49) 



The inequality ( I49> follows from the equation for /?o and the fact that, even for equal leading terms in the sums 
above (with k — K and even «), the difference between Rq and R„ is large due to the next to leading term k — K~ 1, 

Rely(Xf-,))(K-l) 

which is proportional to e t ^ 1. Thus, we arrive at a contradiction with our assumption Rq—R„ 0+, and, 
consequently, it cannot be true. Therefore, for large volumes the real root Ao always gives the main contribution to the 
GCP J42t . and this is the only root that survives in the limit V ^ °°. Thus, we showed that the model with the fixed 
size of the largest fragment has no phase transition because there is a single singularity of the isobaric partition i39\ . 
which exists in thermodynamic limit. 



VIII. FINITE VOLUME ANALOGS OF PHASES 

If K{V) monotonically grows with the volume, the situation is different. In this case for positive value of Re{v) ^ T 
the leading exponent in the rh.s. of (145 > also corresponds to a largest fragment, i.e. to A; = K(y). Therefore, we can 
apply the same arguments which were used above for the case K{V) = K ^ const and derive similarly equations 
(l47>-j48> for /„ and R„. From /„ « ^^fFJ follows that, when V increases, the number of simple poles in J41I I also 
increases and the imaginary part of the closest to the real A-axis poles becomes very small, i.e /„ ^ for n -C K{V), 
and, consequently, the associated decay/formation time T„ w K{V)[nnT]^^ grows with the volume of the system. Due 
to /„ 0, the inequality i49i cannot be established for the poles with n <C K{V). Therefore, in contrast to the previous 
case, for large K{V) the simple poles with n <C K{V) will be infinitesimaUy close to the real axis of the complex 
A -plane. 

From Eq. j48t it follows that 



R„^!^-^ln 
T K{V)b 



(50) 



for \ jj.\^T and K{V) °°. Thus, we proved that for infinite volume the infinite number of simple poles moves toward 
the real A-axis to the vicinity of liquid phase singularity Xi ^ pi{T,jj.)/T of the isobaric partition lIlallTll and generates 
an essential singularity of function ^{V, pi/T) in i40i irrespective to the sign of the chemical potential jj.. In addition, 
as we showed above, the states with Re{v) ^ T become stable because they acquire infinitely large decay/formation 
time T„ in the limit V ^ °°. Therefore, these states should be identified as a liquid phase for finite volumes as well. 

Now it is clear that each curve in Fig. 2 is the finite volume analog of the phase boundary T ~ pL for a given value 
of K{V): below the phase boundary there exists a gaseous phase, but at and above each curve there are states which 
can be identified with a finite volume analog of the mixed phase, and, finally, at Re{v) » T there exists a liquid phase. 
When there is no phase transition, i.e. KiV) ^ K ~ const, the structure of simple poles is similar, but, first, the line 
which separates the gaseous states from the metastable states does not change with the volume, and, second, as shown 
above, the metastable states will never become stable. Therefore, a systematic study of the volume dependence of free 
energy (or pressure for very large V) along with the formation and decay times may be of a crucial importance for 
experimental studies of the nuclear liquid gas phase transition. 

The above results demonstrate that, in contrast to Hill's expectations |6], the finite volume analog of the mixed 
phase does not consist just of two pure phases. The mixed phase for finite volumes consists of a stable gaseous phase 
and the set of metastable states which differ by the free energy. Moreover, the difference between the free energies 



of these states is not surface-like, as Hill assumed in his treatment |6], but volume-like. Furthermore, according to 
Eqs. J44t and ( 145 > . each of these states consists of the same fragments, but with different weights. As seen above for 
the case Re{v) ^ T, some fragments that belong to the states, in which the largest fragment is dominant, may, in 
principle, have negative weights (effective number of degrees of freedom) in the expression for /?„>o ( I44l i. This can 
be understood easily because higher concentrations of large fragments can be achieved at the expense of the smaller 
fragments and is reflected in the corresponding change of the real part of the free energy —R„^qVT. Therefore, the 
actual structure of the mixed phase at finite volumes is more complicated than was expected in earlier works. 

The Hills' ideas were developed further in Ref. |5|, where the authors claimed to establish the one to one correspon- 
dence between the bimodal structure of the partition of measurable quantity B known on average and the properties of 
the Lee-Yang zeros of this partition in the complex g-plane. The starting point of Ref. JsJ is to postulate the partition 
Zg and the probability Pg{B) of the following form 



[dBW{B)e-^-^ ^ Pg{B) ^ ^'^•■^ 



(51) 



where W{B) is the partition sum of the ensemble of fixed values of the observable {B} , and g is the corresponding 
Lagrange multiplier. Then the authors of Ref. |5] assume the existence of two maxima of the probability Pg{B) (= 
bimodality) and discuss their relation to the Lee-Yang zeros of Zg in the complex g-plane. 

The CSMM gives us a unique opportunity to verify the Chomaz and Gulminelli idea on the bimodality behavior 
of Pg{B) using the first principle results. Let us use the equation ( I38> identifying the intensive variable g with A 
and extensive one B with the available volume V' — > V . The evaluation of the rh.s. of ( I38> is very difficult in 
general, but for a special case, when the eigen volume b is small this can be done analytically. Thus, approximating 
.^{B,,X -iri)^ ^(^, A) - /Tj d^{^,X)/dX, we obtain the CSMM analog of the probability (EB 

where we made the Tj integration after applying the approximation for ^(^,A — /tj). Further evaluation of ( 15 2> 
requires to know all possible solutions of the average volume of the system ^^(V) — V — d,^{^^,X)/dX (a = 
{1,2, . . .}). It can be shown |46] that for the gaseous domain v = Re{v) < —2T (see the right panel of Fig. 4) there 
exist a single solution a = 1, whereas for the domain v = Re{v) > 0, which corresponds to a finite volume analog of 
the mixed phase, there are two solutions a = 1,2. In contrast to the expectations of Ref. the probability (15 2t 

p,(y)ir(A,r,M)-E-^4^/[^(«-^)-^i ^ < o, (53) 



1 



has negative derivative for the whole domain of existence of the isobaric partition iF(A , T, ) l46ll . This is true even for 
the domain in which, as we proved, there exists a finite analog of the mixed phase, i.e. for v — Re{v) > 0. Moreover, 

irrespective to the sign of the derivative — the probability i52l cannot be measured experimentally. Above it 

was rigorously proven that for any real £, the IP S'{X,T,jj.) is defined on the real A-axis only for ^{£,,X) — A > 0, 
i.e. on the right hand side of the gaseous singularity Aq: A > Aq. However, as one can see from the equation ( 14 11 1, the 
"experimental" A„ values belong to the other region of the complex A-plane: Re{Xn>o) < 

Thus, it turns out that the suggestion of Ref. |5] to analyze the probability ( 15 1> does not make any sense because, as 
we showed explicitly for the CSMM, it cannot be measured. It seems that the starting point of the Ref. 1 5] approach, i.e. 
the assumption that the left equation ( 15 U gives the most general form of the partition of finite system, is problematic. 
Indeed, comparing ( I50> with the analytical result ( I53l l, we see that for finite systems, in contrast to the major assumption 
of Ref. Isl, the probability W of the CSMM depends not only on the extensive variable V, but also on the intensive 
variable A , which makes unmeasurable the whole construct of Ref. fs*] . Consequently, the conclusions of Ref. \E] on the 
relation between the bimodality and the phase transition existence are not general and have a limited range of validity. 
In addition, the absence of two maxima of the probability (15 3t automatically means the absence of back-banding of 
the equation of state jsj. 



IX. GAS OF BAGS IN FINITE VOLUMES 



Now we will apply the formalism of the preceding sections to the analysis of the Gas of Bags Model (GBM) 1311148 
in finite volumes. In the high and low temperature domains the GBM reduces to two well known and successful 
models: the hadron gas model and the bag model of QGP |50]. Both of these models are surprisingly successful 
in describing the bulk properties of hadron production in high energy nuclear collisions, and, therefore, one may hope 
that their generalization, the GBM, may reflect basic features of the nature in the phase transition region. 

The van der Waals gas consisting of n hadronic species, which are called bags in what follows, has the following 
GCP 1311 
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^2 (^) is the particle density of bags 



where UT)=gk <t>{T,mk) = £ p^dp exp [- {p^ + miyi^/T] 

of mass m^. and eigen volume and degeneracy gi^. This expression differs slightly form the GCP of the simplified 
SMM (|3}, where pL —Q and the eigen volume of ^-nucleon fragment kb is changed to the eigen volume of the bag Vk- 
Therefore, as for the simplified SMM the Laplace transformation with respect to volume of Eq. J54> gives 



Z{sJ) 



-1 



(55) 



7=1 

In preceding sections we showed that as long as the number of bags, n, is finite, the only possible singularities of 
Z{s,T) (I55> are simple poles. However, in the case of an infinite number of bags an essential singularity of Z,{s,T) 
may appear. This property is used the GBM: the sum over different bag states in (I54> can be replaced by the integral, 
Y.J=igj.-- — dmdv...p(m,v), if the bag mass-volume spectrum, p(m,v), which defines the number of bag states in 
the mass-volume region [m,v;m + dm,v + dv], is given. Then, the Laplace transform of Z{V, T) reads llsill 



Zgb{s,T) = / dV e-'"^ Z{yj) ^[s - f{T,s)] ' , where f{T,s) = / dmdv p{m,v) e"''" (j){T,m 



(56) 



Like in the simplified SMM, the pressure of infinite system is again given by the rightmost singularity: p{T) = 
Ts*{T) = T ■ max{sH{T),SQ{T)}. Similarly to the simplified SMM considered in Sect. II and III, the rightmost 
singularity s*{T) of Z{s,T) i56\ can be either the simple pole singularity sh{T) = f {T,sh{T)) of the isobaric 
partition i56\ or the sq{T) singularity of the function f{T,s) ( I56t it-self. 

The major mathematical difference between the simplified SMM and the GBM is that the latter employs the two 
parameters mass-volume spectrum. Thus, the mass-volume spectrum of the GBM consists of the discrete mass-volume 
spectrum of light hadrons and the continuum contribution of heavy resonances jsill 



p{m,v) = ^ gj d{m-mj) d{v-Vj) + &{v-Vo)&{m-Mo-Bv)C v^{m-Bv)^ exp 
i=i 
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(57) 



respectively. Here mj < Mo, vj < Vb, Mq « 2 GeV, Vo « 1 fm^, C, y, 5 and B (the so-called bag constant, B « 400 
MeV/fm^) are the model parameters and 
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is the Stefan-Boltzmann constant counting gluons (spin, color) and (anti-)quarks (spin, color and u, d, i-flavor) degrees 
of freedom. 

Recently the grand canonical ensemble has been heavily criticized 1*521 Is^ . when it is used for the exponential 
mass spectrum. This critique, however, cannot be applied to the mass-volume spectrum (|^} because it grows less fast 
than the Hagedorn mass spectrum discussed in |52, 53] and because in the GBM there is an additional suppression of 
large and heavy bags due to the van der Waals repulsion. Therefore, the spectrum j57t can be safely used in the grand 
canonical ensemble. 

It can be shown ll48il that the spectrum (I57> generates the sq{T) — — j singularity, which reproduces the bag 
model pressure p{T) ~ Tsq{T) fsO] for high temperature phase, and sh{T) singularity, which gives the pressure of 
the hadron gas model |49] for low temperature phase. The transition between them can be of the first order or second 
order or cross-over, depending on the model parameters. 



However, for finite systems the volume of bags and their masses should be finite. The simplest finite volume 
modification of the GBM is to introduce the volume dependent size of the largest bag n — n{V) in the partition J54l l. 
As we discussed earlier such a modification cannot be handled by the traditional Laplace transform technique used in 
is 1 , 48 1 , but this modification can be easily accounted for by the Laplace-Fourier method 1 28 ] . Repeating all the steps 
of the sections V and VI, we shall obtain the equations (|^)-(|^}, in which the function ^(^,A) should be replaced 
by its GBM analog f{l,VB) = ///(A) +/e(A,yB) defined via 

, Vb/Vo 

•Jin r 

fH{X) = Y^ gj^iT,mj)e-'-r\ and /e(A,yB)= V^dk a{T .V^k) e'^'^^-^)^ . (59) 

In evaluating \59\ we used the mass-volume spectrum \51l with the maximal volume of the bag Vg and changed 
integration to a dimensionless variable k = v/Vq- Here the function v) = u{T)v^^^^^ is defined by u{T) = 

The above representation (I59> generates equations for the real and imaginary parts of A„ = /?„ + //„, which are very 
similar to the corresponding expressions of the CSMM (I44> and ( 145 > . Comparing (I59> with (I43> . one sees that their 
main difference is that the sum over k in ( I43> is replaced by the integral over k in j59> . Therefore, the equations (I44> 
and ( 145 1 remain valid for R„ and /„ of the GBM, respectively, if we replace the k sum by the integral for KiV) — Vb /Vq, 
b = Vo, v(A) = Vo{sq{T) — A) and ^k>i{T) = Vo a{T,Vok). Thus, the results and conclusions of our analysis of theR„ 
and /„ properties of the CSMM should be valid for the GBM as well. In particular, for large values of K{V) — Vb /Vq ™d 
Rn < sq{T) one can immediately find out /„ « TTn/Ve and the GBM formation/decay time T„ = yB[7rnrVo]^^ These 
equations show that the metastable A„>o states can become stable in thermodynamic limit, if and only if Vb ^ V . 

The finite volume modification of the GBM equation of state should be used for the quantities which have VAo ~ 1 . 
This may be important for the early stage of the relativistic nuclear collisons when the volume of the system is small, 
or for the systems that have small pressures. The latter can be the case for the pressure of strange or charm hadrons. 



X. HILLS AND DALES MODEL AND THE SOURCE OF SURFACE ENTROPY 

During last forty years the Fisher droplet model (FDM) has been successfully used to analyze the condensation 
of a gaseous phase (droplets or clusters of all sizes) into a liquid. The systems analyzed with the FDM are many and 
varied, but up to now the source of the surface entropy is not absolutely clear. In his original work Fisher postulated 
that the surface free-energy /^i of a cluster of A-constituents consists of surface (A^/^) and logarithmic (InA) parts, i.e. 
Fa = a{T) A^l^ + rrinA. Its surface pait a{T) A^l^ = ao[l - T /Tc\ A^l^ consists of the surface energy, i.e. Oo A^l^, 
and surface entropy —Oq/Tc A^I^. From the study of the combinatorics of lattice gas clusters in two dimensions. Fisher 
postulated the specific temperature dependence of the surface tension (7(r)|FDM which gives naturally an estimate for 
the critical temperature 7^ . Surprisingly Fisher's estimate works for the 3-d Ising model |54], nucleation of real fluids 
|55, 56], percolation clusters |57] and nuclear multifragmentation |3]. 

To understand why the surface entropy has such a form we formulated a statistical model of surface deformations 
of the cluster of A-constituents, the Hills and Dales Model (HDM) i29il . For simplicity we consider cylindrical 
deformations of positive height /zj. > (hills) and negative height — /ij. (dales), with A:-constituents at the base. It 
is assumed that cylindrical deformations of positive height hj^ > (hills) and negative height — /za- (dales), with k- 
constituents at the base, and the top (bottom) of the hill (dale) has the same shape as the surface of the original cluster 
of A-constituents. We also assume that: (i) the statistical weight of deformations exp{—ao\ASk\/si/T) is given by the 
Boltzmann factor due to the change of the surface |A5jt| in units of the surface per constituent si; (ii) all hills of heights 
hk ^ {Hji^ is the maximal height of a hill with a base of ^-constituents) have the same probability dhj^/Hi^ besides 
the statistical one; (iii) assumptions (i) and (ii) are valid for the dales. 

The HDM grand canonical surface partition (GCSP) 
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corresponds to the conserved (on average) volume of the cluster because the probabiUties of hill and dale of the 
same ^-constituent base are identical 12^ 
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Here Pj^ is the perimeter of the cyHnder base. 

The geometrical partition (degeneracy factor) of the HDM or the number of ways to place the center of a given 
deformation on the surface of the A -constituent cluster which is occupied by the set of {nf = 0, 1,2, ...} deformations 
of the /-constituent base we assume to be given in the van der Waals approximation |29]: 
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where si^ is the area occupied by the deformation of fe-constituent base (k = 1,2, ...), Sa is the full surface of the 
cluster, and K,„ax{SA) is the A-dependent size of the maximal allowed base on the cluster. 

The 0(ii^^) -function in {I) ensures that only configurations with positive value of the free surface of cluster are 
taken into account, but makes the analytical evaluation of the GCSP O very difficult. However, we were able to 
solve this GCSP exactly for any surface dependence of K,„ax{SA) using identity (I37> of the Laplace-Fourier transform 
technique Ii28ll : 
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The poles X„ of the isochoric partition are defined by 
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Our analysis shows that Eq. (|5j has exactly one real root Rq = Ao,/m(Ao) = 0, which is the rightmost singularity, i.e. 
Ro > Re{X„yo). As proved in 1 29], the real root dominates completely for clusters with A > 10. 

Also we showed that there is an absolute supremum for the real root Rq, which corresponds to the limit of 

infinitesimally small amplitudes of deformations, Hj^ 0, of large clusters: sup(7?o) = 1 .06009 = Rq = 2 [e*" — l] ' . 
It is remarkable that the last result is, first, model independent because in the limit of vanishing amplitude of 
deformations all model specific parameters vanish; and, second, it is valid for any self-non-intersecting surfaces. 

For large spherical clusters the GCSP becomes Z(5a) ~ 0.38 14 e^ -^^*^*^^"^ ^ , which, combined with the Boltzmann 

'(ToA2/Vr 



factor of the surface energy e 



cluster (y{T) ~ Oq 



1 - 1.06009^ 



, generates the following temperature dependent surface tension of the large 
This result means that the actual critical temperature of the FDM should be 



Tc — C7o/ 1.06009, i.e. 6.009 % smaller in <7o units than Fisher originally supposed. 



XI. STRATEGY OF SUCCESS 

Here we discussed exact analytical solutions of a variety of statistical model which are obtained by a new powerful 
mathematical method, the Laplace-Fourier transform. Using this method we solved the constrained SMM and Gas of 
Bags Model for finite volumes, and found the surface partition of large clusters. Since in the thermodynamic limit 
the CSMM has the nuclear liquid-gas FT and the GEM describes the FT between the hadron gas and QGP, it was 
interesting and important to study them for finite volumes. As we showed, for finite volumes their GCP functions can 
be identically rewritten in terms of the simple poles A„>o of the isobaric partition ( I39> . We proved that the real pole 
Ao exists always and the quantity TAq is the constrained grand canonical pressure of the gaseous phase. The complex 
roots A„>o appear as pairs of complex conjugate solutions of equation ( 143 > . Their most straightforward interpretation 
is as follows: —TRe{X„) has a meaning of the free energy density, whereas bTIm{X„), depending on its sign, gives the 
inverse decay/formation time of such a state. Therefore, the gaseous state is always stable because its decay/formation 
time is infinite and because it has the smallest value of free energy, whereas the complex poles describe the metastable 
states for7?e(A„>o) > and mechanically unstable states for/?e(A„>o) < 0. 

We studied the volume dependence of the simple poles and found a dramatic difference in their behavior for the 
case with phase transition and without it. For the case with phase transition this formalism allows one to define the 
finite volume analogs of phases unambiguously and to establish the finite volume analog of the T — jU phase diagram 
(see Fig. 4). At finite volumes the gaseous phase is described by a simple pole Ao, the mixed phase corresponds to a 
finite number of simple poles (three and more), whereas the liquid is represented by an infinite amount of simple poles 
at IjU I — > oo which describe the states of a highest possible particle density. 

As we showed for the CSMM and GEM, at finite volumes the A„ states of the same partition with given T and jj. 
are not in a true chemical equiUbrium because the interaction between the constituents generates complex values of 



the effective chemical potential. This feature cannot be obtained within the Fisher droplet model due to lack of the 
hard core repulsion between the constituents. We showed that, in contrast to Hill's expectations 01, the mixed phase 
at finite volumes is not just a composition of two states which are the pure phases. As we showed, a finite volume 
analog of the mixed phase is a superposition of three and more collective states, and each of them is characterized by 
its own value of A„, and, consequently, the difference between the free energies of these states is not a surface-like, as 
Hill argued |6|, but volume-like. 

Also the exact analytical formulas gave us a unique opportunity to verify the Chomaz and Gulminelli ideas 01 about 
the connection between bimodality and the phase transition existence for finite volumes. The CSMM exact analytical 
solution not only provided us with a counterexample for which there is no bimodality in case of finite volume phase 
transition, but it gave us an explicit example to illustrate that the probability which, according to Ref. |5] is supposed 
to signal the bimodal behavior of the system, cannot be measured experimentally. 

All this clearly demonstrates that the exactly solvable models are very useful theoretical tools and they open the 
new possibiUties to study the critical phenomena at finite volumes rigorously. The short range perspectives (SRP) of 
this direction of research are evident: 

1. Study the isobaric ensemble and the excluded volume correction for the clusters of the 2- and 3-dimensional 
Ising models, and find out the reliable signals of phase transition on finite lattices. 

2. Widen or refine the CSMM and GMB analytical solutions for more realistic interaction between the constitients. 
In particular, a more realistic Coulomb interaction between nuclear fragments (not the Wigner-Seitz one!) can be 
readily included now into the CSMM and may be studied rigorously without taking thermodynamic limit. 

3. Deepen or extend the CSMM and GMB models to the canonical and microcanonical formulations, and work out 
the reliable signals of the finite system phase transitions for this class of models. 

The major goals for the SRP are (I) to get the reliable experimental signals obtained not with the ad hoc theoretical 
constructs which are very popular nowadays, but directly from the first principles of statistical mechanics; (II) to work 
out a common and useful theoretical language for a few nuclear physics communities. 

However, even the present (very limited!) amount of exact results can be used as a good starting point to build 
up a truly microscopic theory of phase transitions in finite systems. In fact, the exact analytical solution, which we 
found for finite volumes, is one of the key elements that are necessary to create a microscopic kinetics of PTs in 
finite systems. The formulation of such a theory for nuclear physics is demanded by the reality of the experimental 
measurements: both of the phase transitions which are studied in nuclear laboratories, the liquid-gas and hadron gas - 
QGP, are accessible only via the violent nuclear collisions. As a result, in these collisions we are dealing with the PTs 
which occur not only in finite system, but in addition these PTs happen dynamically. It is known that during the course 
of collision the system experiences a complicated evolution from a highly excited (on the ordinary level) state which is 
far from local equilibrium, to the collective expansion of the locally thermalized state and to a (nearly) free-streaming 
stage of corresponding constituents. 

A tremendous complexity of the nuclear collision process makes it extremely difficult for theoretical modeling. 
This is, in part, one of the reasons why, despite a great amount of experimental data collected during last 25 years and 
numerous theoretical attempts, neither the liquid-gas nor the hadron gas - QGP phase transitions are well established 
experimentally and well understood theoretically. It turns out that the major problem of modeUng both of these PTs in 
dynamic situations is the absence of the suitable theoretical apparatus. 

For example, it is widely believed in the Relativistic Heavy Ion community (RHIc) that relativistic hydrodynamics 
is the best theoretical tool to model the PT between QGP and hadron gas because it employs only the equilibrium 
equation of state |58^. Up to now this is just a wishful illusion because besides the incorrect boundary conditions, 
known freeze-out procedure" |59, 60], which are typically used in the actual hydro calculations |58], the employed 
equation of state does not fit into the finite (and sometimes small!) size of the system because it corresponds to an 
infinite system. On the other hand it is known |61J that hydrodynamic description is limited by the weak (small) 
gradients of the hydro variables, which define a characteristic scale not only for collective hydro properties, but also a 
typical volume for the equation of state. 

Above we showed that for finite systems the equation of state inevitably includes the volume dependence of such 
thermodynamic variables as pressure and energy density which are directly involved into hydrodynamic equations. 
This simple fact is not realized yet in the RHIc, but, probably, the chemical non-equilibrium (which is usually 
implemented into equation of state by hand) is, in part, generated by the finite volume corrections of the GCP if 
this is the case, then, according to our analysis of the finite volume GCP functions, it is necessary to insert the complex 
values of the chemical potential into hydro calculations. 



Unfortunately, at present there is no safe recipe on how to include the finite volume equation of state in the 
hydrodynamic description. A partial success of the hybrid hydro+cascade models which might be considered 

as a good alternative to hydrodynamics, is compensated by the fact that none of the existing hydro+cascade models 
was able to resolve the so called HBT puzzles [131 found in the energy range of the Relativistic Heavy Ion Collider. 
Moreover, despite the rigorous derivation 16^ 16511 of the hydro+cascade equations, the hydrodynamic part of this 
approach is suffering the very same problems of the infinite matter equation of state which we discussed above. 
Therefore, further refinements of the hydro+cascade models will not be able to lift up the theoretical apparatus of 
modeling the dynamics of the finite volume PTs to new heights, and we have to search for a more elaborate approach. 

It turns out that the recently derived finite domain kinetic equations j 64^ 65} can provide us with another starting 
point to develop the first principle microscopic theory of the critical phenomena in finite systems. These equations 
generalize the relativistic Boltzmann equation to finite domains and, on one hand, allow one to conjugate two 
(different!) kinetics which exist in two domains separated by the evolving boundary, and, on the other hand, to account 
exactly for the exchange of particles between these domains. (For instance, one can easily imagine the situation when 
on one side of the boundary separating the domains there may exist one phase of the system which interacts with 
the other phase located on the other side of the boundary.) But, first, the finite domain kinetic equations should be 
generalized to the two-particle distribution functions and then they should be adapted to the framework of nuclear 
multifragmentation and the Gas of Bags Model. In doing this, the exact analytical results we discussed will be 
indispensable because they provide us with the equilibrium state of the finite system and tell us to what finite volume 
analog of phases this state belongs. 

Therefore, a future success in building up a microscopic kinetics of PTs in finite systems can be achieved, if we 
combine the exact results obtained for equilibrated finite systems with the rigorous kinetic equations suited for finite 
systems. There is a good chance for the nuclear multifragmentation community to play a very special role in the 
development of such a theory, namely it may act as a perfect and reliable test site to work out and verify the whole 
concept. This is so because besides some theoretical advances and experience in studying the PTs in finite systems, 
the experiments at intermediate energies, compared to the searches for QGP, are easier and cheaper to perform, and 
the PT signals are cleaner and unspoiled by a strong flow. Moreover, once the concept is developed and verified, it can 
be modified and applied to study other PTs in finite systems, including the transitions to/from high temperature QCD 
and dense hadronic matter planned to be studied at CERN LHC and GSl FAIR. Thus, after some readjustment the 
manpower and experimental facilities of nuclear multifragmentation community can be used for a new strategic aim, 
which is at the frontier line of modern physics. 

Such a program, however, requires the coherent efforts of, at least, two strong and competing theoretical groups, 
an access to the collected experimental data and advanced computer facility. Organizationally it will require a very 
close collaboration with experimental groups. Also it turns out that such resources can be provided by the national 
laboratories only. Therefore, we have to find out an appropriate form of national and international cooperation right 
now. Because in a couple of years it will be too late. 
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